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Rayleigh's  Method  Applied  to  a  Conducting  Liquid  Drop  in  the  Presence  of 
a  Point  Charge 

by  Clyde  Morrison 
Richard  Leavitt 


ABSTRACT 


The  method  that  Rayleigh  formulated  in  1879  is  used  to  determine  the 
effect  of  a  point  charge  on  the  natural  frequencies  of  a  charged  drop. 
Explicit  expressions  are  derived  for  the  resonant  frequencies  and  the 
deformation  of  the  charged  drop.  The  results  are  expressed  in  terms  of  sums 
over  particular  types  of  Clebsch-Gordon  coefficients  and  more  complicated 
sums.  Further  work  is  necessary  before  practical  application  of  the  results 
can  be  made;  y>p 
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The  purpose  of  this  work  is  to  extend  our  previous  method  of  solution  to 
problems  involving  drops  in  an  electric  field^1^  to  charged  droplets. 
Consequently,  we  summarize  some  of  our  more  recent  results  on  the  study  of  the 
dynamics  of  liquid  droplets. 

The  methods  employed  here  are  consistent  with  the  method  employed  by  Lord 
( o\ 

Rayleigh  in  his  original  treatment  of  the  dynamics  of  liquid  droplets.  We 
begin  with  a  brief  rederivation  of  Rayleigh's  results  for  uncharged  and 

/  3  \ 

charged  conducting  droplets.  '  This  is  followed  by  a  more  detailed 
derivation  of  the  problem  of  a  charged  droplet  in  the  presence  of  a  point 
charge.  Extensive  attention  is  given  to  this  latter  result,  since  it  appears 
to  be  new.  Many  of  the  results  given  here  may  have  application  to  the 
formation  of  clouds, ^ 4-8 ^  the  physical  aspects  of  drop  formation, ^  the 
stability  of  electrified  surfaces , ^ 10 ^  and  the  electrical  dispersion  of  liquid 
aerosols . ^ 1 1 ^  » 


RAYLEIGH'S  RESULT 

/  2  \ 

In  his  original  paper,  Rayleigh  assumed  that  the  distance,  r,  from  the 
center  to  a  point  on  the  surface  of  the  droplet  can  be  expanded  in  a  Legendre 
series  as 


r(8,t)  *  aQ{t)  +  £  a^(t)  (cos  8)  ,  (1) 

where  we  have  assumed  that  the  drop  is  symmetric  around  the  z  axis  of  the 
drop.  At  this  point  the  z  axis  can  be  chosen  in  any  direction,  but  later, 
when  we  include  the  electromagnetic  energy,  the  electric  field  will  be  assumed 
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along  this  axis.  The  prime  on  the  sum  will  be  used  throughout  to  denote  the 
absence  of  the  k  ■  0  term.  The  volume  of  the  drop  is  given  by 

2ir  ir  r(9) 

V  -  /  d$  /  d0  /  r'2  sin  0  dr’  .  (2) 

0  0  0 

Since  the  drop  is  symmetric  around  z  and  with  u  *  cos0,  we  have 

v=il  j  r3du  #  (3) 

-1 

From  ( 1 )  we  have 


r  »  a. 


‘0  +  3ao  l  \  pk  +  3ao  jJ,  \  ak'  Pk  Pk' 

» 

and  we  use  this  result  in  (3)  to  obtain 

2ir  r,  3 


(4) 


2 

•  *k 


3  ^2a0  +  3  *  2  a0  l  2k+-J  * 


(5) 


where  we  have  used* 


1 

I  V  du  ■ 

-1 


2A+1  U 


• 


(6) 


and  we  have  terminated  the  approximation  at  the  square  of  a^.  We  will  follow 
this  procedure  throughout;  that  is,  we  assume  <<  aQ  and  products  of  more 

than  two  av  are  ignorable  (k  f  0) .  We  assume  the  fluid  is  incompressible  with 

3 

equilibrium  radius  "a"  so  that  the  volume  (5)  is  a  constant,  4xa  /3?  that  is. 


3 

a 


■  *0  *  3  *0 


(7) 


*See  the  appendix  for  a  number  of  useful  relations  involving  Legendre 
polynomials. 
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or 


a  a  a. 


then 


2 

\ 


r  1  r  *  K  1 

t1  +  2  l  2k+J 


*0  k 


2 

\ 


a0  "  4  ’  1/4  £  2k+1 

k 


(8) 


(9) 


The  result  given  by  (9)  serves  as  a  constraint  on  the  variables  a^. 


The  potential  energy,  Ug,  of  the  drop  due  to  the  surface  tension,  y,  is 
the  surface  area  of  the  drop  multiplied  by  the  surface  tension,  or 


2ir 

U  *  y  J '  d$  Jr  sin  8  ds  , 
3  0 


(10) 


where  ds  is  the  arc  length  along  the  surface  given  by 


Then 


.  2  2  2.2 

ds  ■  dr  +  r  d0 


.  r  2  fdr >2 i 1/2 
dS  -  d0  [r  +  [^g]  ] 


(ID 


Using  the  result  (11)  in  (10),  we  have 


U  =  2iry  f  r  [r*  +  (1  -  y2>  (^)2]V2  dy  , 


-1 


dyJ 


(12) 


where  we  have  used  the  relation 


*  -  sin  0  j-  . 
d0  dy 


Expanding  the  integrand  in  (12)  we  have 


Us  -  2ity  [r2  +  -(-~M  }  (§^)2]  dy  , 


(13) 


5 
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0 


y 


and  since  —  does  not  involve  ag,  we  need  consider  no  higher  terms.  The 
integrals  in  (13)  can  be  evaluated  simply  to  give  (see  the  Appendix) 


U  .  2wy  [2a?  .  I'  a2 

s  L  0  £  2k+1  k 

and  if  we  use  the  constraint  given  in  (9),  we  have 


<]  • 


r-  2  .  r'  (k-1)  (k+2)  2i 

°a  *  2’*  t2a  +  I  — »7i - ' 


which  gives  the  potential  energy  due  to  surface  tension  correct  through  terms 
of  order  a^. 

To  calculate  the  kinetic  energy,  T,  we  need  to  evaluate  the  integral 

T  *=  /  V2  P  v2  dr  ,  (16) 

where  p  is  the  density  and  v  is  the  velocity.  Since  we  are  assuming  that  the 
fluid  is  incompressible  and  that  there  are  no  sources  or  sinks  within  the 
drop,  we  have 


and  further,  if  we  assume 


7  .  v  »  0  ; 


V  x  v  *  0  , 


then  v  can  be  derived  from  a  potential  function,  such  that 

v  **  7ip  , 
v 


and  from  (17)  we  have 


Using  (19)  in  (16)  we  have 


7tJ>  =  0 
v 


T  •  V2  P  /  (Vi|»  )  dx  , 


6 
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and  by  converting  the  integral  over  the  volume  to  a  surface  integral,  we  have 


T  »  y/2  P  /  (V^)  *  <*o 


(22) 


Assuming  that  the  area  element  do  is  approximately  along  r  (corrections  are  of 

♦  A  2 

higher  order),  we  can  write  do  *  r  a  d$  du,  where  <p  «  the  azimuthal  angle 
and  r  *»  the  unit  vector  along  the  radius;  then 


T  *  irpa^  / 
-1 


♦ 

T 


dij> 

Tv 

dr 


dy 


(24) 


The  solution  of  (19)  appropriate  for  our  problem  is 


*  £  8  r  P  (y)  ,  (25) 

rv  *•  n  n 
n 

where  we  have  assumed  that  the  potential  is  evaluated  at  the  surface.  Using 
(25)  in  (24)  we  have 


T 


n  2n-1 
2n+1  a 


(26) 


The  ,8  in  (26)  can  be  evaluated  by  equating  the  velocity  at  the  surface,  given 

• 

by  — —  ,  to  r  from  ( 1 ) ,  or 
or 


Then  (26)  becomes 


•  n-1 

a  »  n  a  8  . 

n  n 


T 


•2 

a 


n 

n( 2n+l ) 


(27) 


(28) 


The  results  given  in  (28)  and  in  (15)  are  sufficient  to  form  the  Lagrangian 

(21 

(L  »  T  -  U)  in  the  absence  of  other  energy  sources.  Rayleigh'  '  uses  these 


7 
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results  to  obtain  the  equation  of  motion  for  the  «n(t) *by  treating  the  an  as 
generalized  coordinates.  We  will  put  off  doing  this  until  the  electromagnetic 
energy  is  contained  in  the  Lagrangian. 


k-  . 


The  particular  problem  we  wish  to  solve  here  is  for  a  charged  conducting 
sphere  of  equilibrium  radius  a  with  total  charge  Q.  This  problem  was  also 
solved  by  Rayleigh, * 12  *  but  by  a  slightly  different  technique  than  was  used  in 
his  original  paper.  Since  we  are  interested  here  in  extending  his  original 
technique  to  a  new  problem,  we  shall  include  more  detail  than  in  our  previous 
discussion. 

From  Maxwell's  equations  we  have 

7  X  E  «  0  ,  (29) 


or  e  *  -  7  if>  . 

Then  from 

7  *  E  -  0  ,  (30) 

we  have 

7^  ip  «  0.  (31 ) 

Thus,  the  potential  is  a  solution  of  Laplace's  equation.  The  appropriate 
solution  of  (31)  for  our  problem  for  r  >  r(0,t)  is 


if  -  I 

n 


A 

— —  p 

n+1  n 
r 


(32) 


The  total  charge  on  the  sphere,  Q,  is  assumed  constant;  thus,  Aq  =  Q.  This 
result  (Aq  -  Q)  is  simply  proven  by  assuming  a  very  large  sphere  surrounding 
the  charged  drop  and  employing  Gauss'  theorem  at  the  surface  of  that  sphere. 

The  other  A  must  be  found  from  the  boundary  conditions  of  the  problem.  Once 

°  (13) 

the  A  have  been  determined,  the  electromagnetic  energy  is  given  by' 
n 


4 
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UE  “  V2  f  p  da  r  ( 33 ) 

where  p  is  the  surface  charge  density,  if>  is  the  potential  at  the  surface,  and 
the  integral  covers  the  surface  of  the  sphere.  Since  the  potential  on  a 
conducting  surface  is  constant,  V,  say,  then  (33)  becomes 

Ue»V2QV  /  (34) 

where  Q  =  /  pda  and  is  a  constant  by  assumption.  The  potential  at  the 
surface,  V,  is  however  a  function  of  the  variation  of  the  shape  of  the 
surface;  as  will  become  apparent. 


and 


To  proceed  further  it  is  convenient  to  introduce  the  notation 

r(0,t)  *  aQ  +  S  l  ak  , 

A  -  [  8v  Alu)  . 
n  *•  n 


(35) 


which  serves  to  keep  track  of  the  order  of  the  corrections .  In  the  final 
result  we  let  6  »  1 ,  since 

6  only  serves  as  an  artifice  to  keep  the  terms  in  order.  The  single  boundary 
condition  needed  is 


«  .  A 
i ^  ! S-f  . 

r  u  n+1  n 


(36) 


n  r 


(v) 


where  r  is  the  function  of  9  given  in  (35).  The  coefficients  A  in  (35)  can 

n 

be  determined  by  using  (35)  in  (36)  and  equating  powers  of  &  on  each  side.  If 

2 

this  is  done,  the  results  through  5  are 


.(0) 


9 


h* 

L:  •. 


8 

fe'-V 

t\'r 

k' 


(1)  „  n-1 

n  3  2  a0  an  ' 
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u 


■‘-t 


f 


v(1)  -o  , 


and 


or 


v(2)  =  (n+1 )  anAn 


(1) 


n+2  < 2n+1 


ao  n 


-  +  fi-  y*  Ja 

)  3  1  2n+1 


.(2) 


_  n  a 

Q  r  _ n 

“  3  \  ( 2n+1 ) 


(37) 


ao  n 


To  obtain  this  result,  the  equations  given  in  the  appendix  were  used.  The 
electromagnetic  energy  through  second  order  is  given  by  using  (37)  in  (34)  or 


°P  »V2Q  [v<0)  +  V(2)]  . 


(38) 


,(0) 


The  result  for  Vv  '  can  be  reduced  further  by  using  the  constraint,  (9),  in 
(37)  to  give 

V(0)  --2  [l  +i-  l'  — .  (39) 

a  k  9  *•  9r>4-1  J 


2 

a  n 


The  result  (39)  along  with  (15)  and  (28)  is  sufficient  to  form  the  Lagrangian 
for  the  system,  that  is 


resulting  in 


L  *  2irpa3  y 


•2 

a 

n 


n( 2n+1 ) 


-  2ry  l 


n 


'  (n-1 ) (n+2)  2 

(2n+1)  an 


(40) 


♦  -V  Z  ' 


(n-1)  2 

.  3  L  2n+1  n 
2a  n 
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where  the  constant  terms  in  (15)  and  in  (39)  have  been  ignored. 


The  equation  of  motion  for  an 
which  becomes 


is  Lagrange's  equation, 


d  3L 


where 


2 


(it 

n 


a 

n 


+  u 


a 


n 


0  , 


atfizli  [Y(n+2)  .  -sL] 

paJ  4itaJ 


(41) 


/  1  \ 

which  is  the  result  derived  by  Rayleigh1  *•'  expressed  in  the  same  form.  As 

2 

was  pointed  out  by  Rayleigh,  the  system  becomes  unstable  for  <  0  ,  which 
occurs  for  values  of  Q  such  that 


>  y(n  +  2)  ,  (42) 

4ira 

2 

which  for  the  lowest  mode  (n  *  2)  requires  that  -*■-  >  4y.  The  charge  on  the 
drop  acts  so  as  to  reduce  the  surface  tension  so^^hat  only  the  higher  modes 
are  stable.  The  result  given  in  (42)  has  been  used  by  numerous  workers  as  a 
starting  point  for  the  discussion  of  not  only  charged  drops  but  also  uncharged 
drops  in  an  external  uniform  electric  field,  one  of  the  more  recent  being  the 
article  by  Smith. 


If  in  (42)  we  replace  the  charge  by  the  electric  field  at  the  surface  of 
the  equilibrium  sphere  (Q  =  a  E),  then  we  can  write 

E2a 

>  n  +  2  ,  (43) 

and  for  the  fundamental  mode  (n  *  2)  we  have 


f 


11 
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for  the  onset  of  instability.  Frequently,  in  reporting  experimental  results 
of  drops  in  electric  fields,  the  field  E  in  (43)  is  taken  as  the  external 
field  and  the  constant  Cg,  determined  in  the  equation 

E  ^  =  CQ  ,  (44) 

by  experimental  means.  Such  a  procedure  was  adapted  by  Wilson  and  Taylor  *' 
in  their  experiments  with  soap  bubbles,  where  they  evaluated  Cg  =  1  .61  ± 
as  the  critical  value  in  (44). 


CONDUCTING  DROP  IN  THE  FIELD  OF  A  POINT  CHARGE 

We  assume  as  before  that  the  conducting  drop  of  charge,  Q,  equilibrium 
radius  a  is  at  the  origin  of  the  coordinate  system.  The  point  charge,  Q1 ,  is 
located  on  the  positive  z  axis  at  a  distance  R  from  the  origin.  The  energy  of 
the  system  corresponding  to  (33)  is 

UE  "  2  27  +  2  Qi*(R)  '  (45) 

where,  as  before,  V  is’ the  potential  at  the  surface  of  the  sphere  and  i|»(R)  is 
the  electric  potential  at  the  point  charge.  The  expressions  for  the  kinetic 
energy,  (28),  and  the  surface  energy,  (15),  remain  the  same  as  before. 

The  potential  ip  appropriate  for  the  problem  is 
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f" 

! 


fa 


where  in  the  second  term  we  have  assumed  r  <  R,  which  is  the  region  of  the 
interest.  It  is  convenient  to  rewrite  (46)  as 


•  Q 


*  -  £  +  +  J'  -2-  P  +  Q  V’  P  , 

y  r  R  n+1  n  1  L  _n+1  n 


(47) 


n  r 


so  that  all  the  sums  have  the  terms  n  *  0  omitted.  To  evaluate  the  constants 
in  (47),  it  is  convenient  to  assume  the  same  "bookkeeping  procedure"  as  given 
in  (35)  to  expand  each  side  of  the  equation 


■  *£  +  iT+  +  81  l 


(48) 


and  equate  the  coefficients  of  each  order.  The  terms  in  the  potential  are 
then  given  by 


.(0) 


n-1 
n  a„  a 


♦  f-  , 


(n+1 )a  A(0) 
n  n 


( 2n+1  )R 
(1  ) 


(2n+1  )  a, 


n+2 


n  a  (n+1 )a  A 

(2)  _  _2_  y  n  r  n  n 


(n+1 ) (n+2)A 


(0) 


(2n+1 )a 


n+2 


n  2(2n+1 )a 


n+3 


n-2 


+  Q  T  Ty- — 7-r - —  7  a  a.  ,<kk'  n>  , 

*1  L  2(2n+1)  _n+1  ,  k  k'  1 

n  R  kk ' 


(49) 

(50) 


l  <kk'|n>' 

k,k' 


(51 ) 


where  the  Clebsch-Gordan  coefficients  (C-G  coefficients)  <kk’|n>'i  are  given  in 
the  appendix.  The  corresponding  expression  for  the  An  are  given  by 

Q, 


(0) 


‘1  2n+1 

- T  a„ 

.n+1  0 


(52) 


(1) 


Aa 


,  ,  ’  U+1)A(o0)a. 

A'"  -  a”  Q  a  +  a"”  \  - ; -  <kijn>  -  a"  Q  \ 

n  0  n  0  k  A  0  1  *k  RU1 


13 


°  a^  <kJt|n>2 
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(53) 


(2) 


_  ,  ,  U+DA^a. 

-ao  8£AaK'<  ki',|n>  +  *o  l— - <ktin> 

*Jc  a_ 


kk1 


-  a 


(£+1 )  (£+2)a[°  A(A-1)a^Q, 

y  r - L_  + _ 2_L 

o  .L,  L  0_a+i  _ Jt+i 


n-1 


£kk' 


2a 


2R 


]  *kak,  l  <kk' | v>2<vi|n>2  . 


By  combining  (49)  through  (54),  we  obtain 


.(0) 


11  + 

R  ao  ' 


.( 1 )  r  r  O^n+I 


4 1  (irY 


.(2) 


n 

*0  2  ‘  . 

n  a  Q  a.  n+1  _ 

3  ^  2n+1  +  T  Vifc")  <kA  )n>  , 


ao  n 


(1) 


n-1 


n 
*  ( 2) 


a.  nki 

<2t-H)aJ+1n 

n-1  .  r  0_i<la|n>2  _ 


-  ao  '  A  -  ao  '  a,  I  —J?, 

ik  R 


(55) 

(56) 

(57) 

(58) 


r  ■  i  **,%  «!->*  -  vr  tL.!2t+,,t^+,  vk-v*  • 


Ak 


and 


Tk'ki  *  £  v  <k'£lv>  <vk|n>  .  (59) 

v 

These  results  are  all  that  are  necessary  to  obtain  the  electrostatic  energy 

2 

through  terms  of  order  ak«  In  (56)  through  (58),  Sq  can  be  replaced  by  "a", 
but  in  (55)  the  constraining  equation,  (9),  must  be  used  to  cast  this  term  in 
final  form.  If  this  is  done  we  have  V  ■  ^  and  if  W1  =  1/2 


QV,  the  first  part  of  the  energy  given  in  (45)  is 


1 
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I 


w 


2  p.  +  2  +  II  i 

5  Lr  a  2  ^ 


n+1 


1  2  LR 


a  2 

a  n 


_  _J>  £  (n-1  )  2 


3  2n+1  -  .  . 

an  a  ki 


‘n  tat\Gkl(*)J 


where 


1 


and 


G^U)  »  l  xn+1  <  ki|n>' 
*"  n»1 

x  ■  a/R 


N 


M 


This  result,  (60),  is  in  final  useable  form  for  Lagrange’s  equations,  but  the 
second  part  of  the  electrostatic  energy  given  by  (45)  needs  more  reduction 
before  it  can  be  used.  From  the  second  term  in  (45)  we  have 


»2  -  i «,  *<»>  • 


(61) 


where  the  terms  in  i|»(r)  that  become  infinite  at  r  =  R  arc  removed.  From  (46) 
and  (48)  we  have 


(62) 


n  R 


where  in  the  second  term  we  have  used  the 

relation  P  (cosQ)  *  1  for  0*0,  since  we  are  assuming  that  the  charge  Q.,  lies 
n  * 

on  the  positive  z  axis  at  a  distance  R  from  the  origin.  As  in  our  previous 
procedure,  we  assume 


W2  -  w‘0)  *  w'”  t  W<2)  . 


(63) 


so  that  from  (62)  we  have 


(0) 


.(0) 


[$♦  I'-fcrJ  . 


(64) 


n  R 
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,(D  1 


1  .  r  n 

2  21  ^  _n 


n  R ' 


,( 2)  1 


1  «  c'  n 

2  C1  ^  _n+1 


Each  of  the  terms,  (64),  (65),  and  (66),  require  considerable  effort  to  reduce 
them  to  useable  results. 


The  first,  (64),  is  the  simplest;  when  (52)  is  used,  we  have 


w(0>  „  1  ?1.  °1  r  x^(3-x“)  1  r 1  “n  ] 

W2  2  R  2a  L  2  2,2  2  *  2n+1  J  '  (67) 

1 -x  ( 1  -x  )  a  n 

where  x  =  a/R,  and  we  have  used  the  constraint  (9).  If  in  (65)  we  use  (58), 
we  have 


x4<3-x2)  1 


"2”  -J?  l' 


I  *1  n 1  lr+1 

a - -  l  (2k+1 )  x*  a  G.  (x)  ,  (68) 


2a  n,k 


where  G^fx)  is  as  defined  in  (60).  The  last  term,  (66),  can  be  obtained  by 
substituting  (59)  into  (66),  which  gives 

(2)  ®®i  r'  . 

W2  -—I  Aa,ak  Gffc(x) - 7  l  m(2m+1)akal  Gfcm(x)  G^tx)  .(69) 


2a3  i 


3  L 
2a  k£m 


With  the  result  given  in  (69),  we  have  completed  all  the  terms  in  the 
potential  energy  that  are  necessary  to  form  the  Lagrangian  of  the  problem. 
Unlike  most  of  the  previous  results,  there  appears  to  be  no  advantage  in 
combining  (67),  (68)  and  (69),  since  corresponding  terms  do  not  exist.  The 
Langrangian  can  be  written 

T  -  U  ,  (70) 
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(71) 


where  T  is  given  in  (28),  Ug  is  given  in  (15)  and 


u  -  w,  *  w'°»  +  w'1 >  ♦  w<2> 
E  1  2  2  2 


(72) 


The  equation  of  motion  for  the  generalized  coordinate  a  is  given  by 


d  3L 

dt  • 
3a 


3L 

3a 


(73) 


which  can  be  written  as 


a _ jw 

dt  +  3aP  ' 
With  the  appropriate  substitutions  we  have 


30. 

3a 


E 


(74) 


4  3  - E -  +  4~y  lezilLEtj)  . 

p  p(  2p+1  )  +  4lrY  2p+1  ap 

'  a2  (1-x2)2  2P*’ 

+  ~ r  a  Tm  (2m  +  1)  G2  (x) 

a3  p  m  P" 


+  P0  +  F1  ' 


(75) 


where 
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QQ.  ryi.1  2.  _l  V+1 

— 2  +  — 2  E  (2k  +  1)  x  Gk  (x)  » 

ar  2a  k  p 


««  p  • 

- —l  “  *  p)  at  V 


U1  P*  P* 

4  — —  T  T  m(2m  +  1)  a.  G.  (x)  G  (x) 

ILL  l  tm  nn 


a  l  m 


The  terms  on  the  right  side  of  (75)  have  been  grouped  in  a  particular 
manner.  The  first  four  terms  multiply  ap  and  act  as  a  restoring  (or 
repulsive,  depending  on  the  sign  of  )  force.  The  term  Fq  acts  as  a  constant 
force  causing  a  displacement  from  the  equilibrium  point  (normally,  the 
equilibrium  point  is  at  ap  =  0),  and  the  term  F1  acts  to  couple  ap  to  the 
other  modes  of  vibration. 


If  we  ignore  FQ  and  F1  on  the  right  side  of  (75),  we  can  write 


a  +  a)  a  *  0  , 

P  P  P 


where 


4  2. 

pa3a>2  *  Yp(p  -  1)(p  +  2)  -  U  p(p  -  1)  +  U  - - T  P 

p  1  2  (1-xV 

0  2 
+  0,„(2p  )(2p  +  1)  G  (x)  -  0,  p(2p  +1)  l  m(2m  +  1)G  (x)  , 

12  pp  2  "  Pm 


r  °2 


T  and  °J2 


The  first  two  terms  on  the 


o  1 

where  0  *  *  — ,  U  *  - r,  and  U  ,  ®  - -  .  The  first  two  terms  on  the 

right  side  of  (77)  area  just  pa3  times  the  frequency  of  a  charged  drop  as 
derived  originally  by  Rayleigh, ^ 1 and  are  given  in  (41).  Since  Fq,  F., ,  U2» 
and  0^2  vanish  when  ■  0  the  result  (77)  becomes  identical  to  (41),  as  it 
should.  If,  on  the  other  hand  Q  -  0,  we  have,  from  (77), 
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which  can  cause  either  an  increase  or  decrease  in  the  resonant  frequency  of  a 
3  2 

free  drop  [pa  u  =  YP(P  -  1)(p  +  2)1  through  the  rather  complicated 

P 

dependence  of  the  last  two  terms  on  x.  The  frequency  shift  given  in  (78)  does 
not  depend  on  the  sign  of  the  point  charge  Q1  as  it  should  not;  however,  when 
the  charge  Q  is  not  zero,  then  (77)  shows  that  the  sign  of  the  charge  Q1 
relative  to  Q  is  significant,  as  apparent  by  the  terms  involving 

Many  of  the  less  obvious  results  of  (77)  will  have  to  await  considerable 
computational  investigation  before  quantitative  statements  can  be  made. 
However,  it  is  possible  that  further  analytical  work  can  be  done  on  the  Gpp(x) 
and  the  sums  involving  these  functions.  They  seem  to  be  expressible  in  terms 
of  terminating  hypergeometric  functions.  This  can  be  seen  by  expressing  the 
CG-coeff icients  in  Gpp^x)  in  Racah's^®'1^  closed  form  and  using  expressions 
for  the  hypergeometric  series  and  their  identities  given  in  Rainville.^ 17 ^ 


APPENDIX 


During  the  course  of  the  development  of  the  theory  a  number  involving 
relations  of  Legendre  polynomials  are  needed.  We  abbreviate  P  ( y )  =  P  and 

^  X  *  X 

—  P^(y)  **  P£  in  the  results. 


-1 


24+1  4k 


r  2.  ,  24(4+1)  , 

/  (1  -  v  )  P'.  Pv  <*y  =  - - -  <$ 


-1 


/  (1  -  U  )  PnP'P^  du  = 


-1 


4  k  H  24+1  Ik 

[4(4+1  )+k(k+1  )-n(n+1  )  ] 
2n+1 


<4k  n>‘ 


(A1) 

(A2) 

(A3) 


[<4k|n>  =  <4(0)k(0) | n{0)>  *  C  (4kn;00)]  , 


where  <4k|n>  is  a  Clebsch-Gordan  coefficient.  Of  these  three  results,  the 
last  (A3)  is  the  only  difficult  one  to  derive.  A  thorough  discussion  of 
Clebsch-Gordon  coefficients  is  given  by  Rose^1*^  or  Brink  and  Satchler, ^ 1 
but  is  more  general  than  is  necessary  here.  In  general,  the  simple  relations 
given  here  can  be  derived  from 


Vl  *  l  <ik|n>2  Pn 


(A4) 


Thus, 


<4k  n> 


(2n+1) 

2 


/  P  P.P,  dy 
J  n  k  l 


-1 


Mg) 

which  is  a  special  case  of  Gaunt's'  '  formula. 


Other  properties  which  we  need  are 


(A5) 


<4k|n>^  «  <kl|n>^  , 


(A6) 
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<lk  n> 


<41  n> 


<lk|o>' 

1+1 

21+1  ‘ 

_i _ 

21+1  ' 


0  . 


I>2n+1 

^2k+1 


)  <ln|k> 


2 


t 


a  -i—  —  5 

21+1  Ik 


n  -  1  +  1 
n  -  1  -  1 


For  other  values  of  n 


(A7) 

(A8) 


(A9) 


The  result  given 'in  (A4)  can  be  used  to  reduce  more  complicated  products  such 
as 


P^P.P  *  J  <kl|n'>  P  ,P 
kin  **  1  n  n 

n’ 

l  <kl|n‘>2<n,n|n">2PnB  , 


(A10) 


n\n" 


which  gives  the  result  (A5)  if  both  sides  are  integrated  over  y,  and  (A1 )  with 
(A8)  is  used. 
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